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We calculate the partition function Z(G, Q, v) of the Q-state Potts model exactly for strips of the 
square and triangular lattices of various widths L y and arbitrarily great lengths L x , with a variety 
of boundary conditions, and with Q and v restricted to satisfy conditions corresponding to the 
ferromagnetic phase transition on the associated two-dimensional lattices. From these calculations, 
in the limit L x — > oo, we determine the continuous accumulation loci B of the partition function 
zeros in the v and Q planes. Strips of the honeycomb lattice are also considered. We discuss some 
general features of these loci. 
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I. INTRODUCTION 

This is the second in a series of two papers on zeros of 
the Q-state Potts model partition function Q-Q on lat- 
tice strip graphs of fixed width L y and arbitrarily great 
length L x , with Q and the temperature- like variable v 
restricted to satisfy the condition for the ferromagnetic 
phase transition on the associated two-dimensional lat- 
tice. From these calculations, in the limit L x — > oo, we 
exactly determine the continuous accumulation loci B of 
the partition function zeros in the v and Q planes. These 
loci are determined by the equality in magnitude of the 
eigenvalues of the transfer matrix of the model with maxi- 
mal modulus and hence are also called the set of equimod- 
ular curves (where "curve" is used in a general sense that 
also includes line segments). In the first paper Q we 
carried out this study for self-dual strips of the square 
lattice and found a systematic pattern of features, from 
which we were able to conjecture properties applicable 
for strips of arbitrarily large widths. Here we continue 
this study by considering strips of the square, triangu- 
lar, and honeycomb lattices with a variety of different 
boundary conditions and studying how these boundary 
conditions affect the loci B. 

We begin by briefly recalling the definition of the model 
and some relevant notation. On a graph G at tempera- 
ture T, the Potts model is defined by the partition func- 
tion 

Z(G,Q,v) = ^2e^p(K^2S aiaj ) (1.1) 



where cr^ = 1 , . . . , Q are the classical spin variables on each 
vertex (site) i € G, (ij) denotes pairs of adjacent vertices, 
K = (3J where (3 = (fc^T) -1 , and J is the spin-spin 
coupling. We define a — e K and v = e K — 1, so that v has 
the physical range of values < v < oo and — 1 < v < 
for the respective ferromagnetic and antiferromagnetic 



cases J > and J < 0. The graph G = G(V,E) is 
defined by its vertex set V and its edge (bond) set E. The 
number of vertices of G is denoted as n — n(G) = \V\ 
and the number of edges of G as e(G) = \E\. The Potts 
model can be generalized from non-negative integer Q 
and physical v to real and, indeed, complex Q and v via 
the cluster relation @ Z(G, Q,v) = Eccg Q fc(G V( G '), 
where G' = {V, E') with E' C E, and fc(G') denotes the 
number of connected components of G' . 

The Potts model partition function is equivalent to an 
important function in mathematical graph theory, the 
Tutte polynomial T(G,x,y) [IHH: 

Z(G, Q, v) = (x- l) k ^(y - 1)"( G > T(G, x, y) (1.2) 



where 



Q 



(1.3) 



The phase transition temperatures of the ferromag- 
netic Potts model (in the thermodynamic limit) on the 
square (sq), triangular (t), and honeycomb (he) lattices 
are given, respectively, by the physical solutions to the 
equations P,P 



Q 



(sq) 



and 



Q = v 2 {v + 3) (i) 



3Qv - v 3 = (he) 



(1.4) 



(1.5) 



(1.6) 



The conditions ljl.5|l and ljl.6|l are equivalent, owing to 
the duality of the triangular and honeycomb lattices. In 
terms of the Tutte variables, these conditions are (y — 



l)(x 
A = t. 



y) = for A = sq, (y 
and (y — l) 2 (x 2 



i)(z/ 2 
y-i) = 



\-y-x- 
for A 



1) = 

= he. 



-- for 
Since 



y = 1 means K — 0, i.e., infinite temperature, the root at 
y = 1 is not relevant; dividing both sides of these three 
equations by the appropriate powers of (y — 1), we thus 
obtain the conditions x = y (sq), y 2 + y — x — 1 = 
(t), and x 2 + x - y - 1 = [Hj. 

Although the infinite-length, finite-width strips that 
we consider here are quasi-one-dimensional systems and 
the free energy is analytic for all nonzero temperatures, it 
is nevertheless of interest to investigate the properties of 
the Potts model with the variables Q and v restricted to 
satisfy the above conditions. For a L y x L x section of the 
respective type of lattice, as L x — > oo and L y — > oo with 
Ly/L x equal to a finite nonzero constant, one sees the 
onset of two-dimensional critical behavior. The strips 
with L x — ► oo and L y fixed provide a type of interpo- 
lation between the one-dimensional line and the usual 
two-dimensional thermodynamic limit as L y increases. 
One of the interesting aspects of this interpolation, and 
a motivation for our present study, is that one can ob- 
tain exact results for the partition function Z(G,Q,v) 
and (reduced) free energy / = liru n ->oo n ~ 1 m The 
value of such exact results is clear since it has not so 
far been possible to solve exactly for f(A,Q,v) for arbi- 
trary Q and v on a lattice A with dimensionality d > 2, 
and the only exact solution for arbitrary v is for the 
d = 2 Ising case Q — 2. Thus, exact results on the 
model for infinite-length, finite-width strips complement 
the standard set of approximate methods that are used 
for d > 2, such as series expansions and Monte Carlo 
simulations. Although the singular locus B v does not in- 
tersect the real axis on the physical finite-temperature 
interval — 1 < v < oo for the infinite-length, finite- width 
strips under consideration here, properties of the image 
of this locus under the above mappings H1.4 |) -l|1.6 [) . Bq, 
can give insight into the corresponding locus Bq for the 
respective two-dimensional lattices. Indeed, one of the 
interesting results of the present work and our Ref. Q 
on infinite-length, finite- width strips is the key role of the 
point Q = 4 for the locus Bq , which can make a connec- 
tion with the locus Bq for the physical phase transition 
of the Potts model on two-dimensional lattices. 

We next describe the boundary conditions that we con- 
sider. The longitudinal and transverse directions of the 
lattice strip are taken to be horizontal (in the x direction) 
and vertical (in the y direction). Boundary conditions 
that are free, periodic, and periodic with reversed orien- 
tation are labelled F, P, and TP. We consider strips 
with the following types of boundary conditions: 

1. (FBC y , FBC X ) = free 

2. (FBC y ,PBC x ) = cyclic (eye.) 

3. (FBC y ,TPBC x ) = Mobius (Mb.) 

4. (PBC y ,FBC x ) = cylindrical (cyl.) 

5. (PBC y ,PBC x ) = toroidal (tor.) 

6. (PBC y ,TPBC x ) = Klein-bottle (Kb.) 



We thus denote a strip graph of a given type of lat- 
tice A = sq or t as A[L y x L X ],BC, where BC =frce 
for (FBC y , FBC X ) and similarly for the other boundary 
conditions. In earlier work we showed that although the 
partition functions Z(A[L y x L x ], BC,Q,v) are different 
for cyclic and Mobius boundary conditions, B is the same 
for these two, and separately that although this partition 
function is different for toroidal and Klein-bottle bound- 
ary conditions, B is the same for these two latter condi- 
tions. Therefore, we shall focus here on the cases of free, 
cyclic, cylindrical, and toroidal boundary conditions. 

Our procedure for calculating Z(G,Q,v) on these 
strips is as follows. For the square and triangular lat- 
tices, the equations <|1.4I) and 11.5fl have the simplifying 
feature that they have the form Q = gh_{v), where g^(v) 
is a polynomial in v. Accordingly, to restrict Q and v to 
satisfy the phase transition conditions for the respective 
two-dimensional lattices, we start with the exact parti- 
tion function and replace Q by ghiv) for A = sq and 
A = t. We then solve for the zeros of Z(G, gA{v),v) and, 
in the L x — > oo limit, the continuous accumulation loci 
B v in the complex v plane. The image of these zeros and 
loci under the respective mappings (|1.4f> and l|1.5f) yield 
the zeros and loci in the complex Q plane. In the case 
of the honeycomb (he) lattice, the ferromagnetic phase 
transition condition (|1.6|) is nonlinear in both v and Q. 
Since it is of lower degree in Q, we solve for this vari- 
able, obtaining Q — (v/2)(— 3 ± y9 + Av ). Of course, 
only one of these solutions is physical for the actual two- 
dimensional lattice, namely the one with the plus sign. 
Given the fact that the triangular and honeycomb lat- 
tices are dual to each other and the consequence that 
properties of the phase transition of the ferromagnetic 
Potts model on the triangular lattice are simply related 
to those on the honeycomb lattice, it follows that, inso- 
far as we are interested in applying our exact results on 
infinite-length, finite-width strips to two-dimensions, it 
suffices to concentrate on strips of either the triangular 
or honeycomb lattice. Since the condition (|1.5|l is easier 
to implement than the solution for Q given above, we 
shall mainly focus on strips of the triangular lattice, but 
also include some comments on honeycomb-lattice strips. 

We denote the Tutte-Beraha numbers 0, 0, EI 



Q r = 4cos 2 (7r/r) . (1.7) 



For the range of interest here, 1 < r < oo, we note that 
Q r monotonically decreases from 4 to as r increases 
from 1 to 2, and then Q r increases monotonically from 
to 4 as r increases from 2 to oo. For our analysis of 
strips of the square lattice, it will also be useful to denote 
v r = — 2 cos(7r/r) so that Q r — v 2 .. Further background 
is given in Ref. |^. 
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II. SOME GENERAL STRUCTURAL 
PROPERTIES 

For these strips, the partition function has the general 
form (with L x = m) 

Z(A[L y x m],BC, Q,v) = £ Cj (X A ,BC,L v>j ) m (2.1) 



where the coefficients Cj are independent of m. It will be 
convenient to separate out a power of v and write 



^A,BC,L y ,j — v Ly \a,bc,l v 



(2.2) 



We denote cyclic strips of the lattice A of width L y 
and length L x = m as by A\L „ x m], eye.. The partition 
function has the general form [13, Il3l [lj] 

L y n z (L y ,d) 
d=0 3=1 

(2.3) 

where we use simplified notation by setting 

^A,cyc.,L y ,d,j = ^-A,L y ,d,j and \A,cyc.,L y ,d,j = ^A,L y ,d,ji 

and where 



nz(L y , d) 



{2d+l) f 2L y 

(Ly+d+1) \Ly-d 



for < d < L y and zero otherwise, and 



d 
j=0 



2d- j 



Q 



d-3 



(2.4) 



(2.5) 



The coefficients can also be expressed in terms of 
Chebyshev polynomials of the second kind as c&> = 
U2d{^-)- The first few of these coefficients are c' ' = 1, 
C W =Q-1, = Q 2 - 3Q + 1, and = Q 3 - 5Q 2 + 
6Q — 1. The form H2.3J) applies for cyclic strips of not 
just the square lattice, but also the triangular and hon- 
eycomb lattices [lj, [l5j . The total number of eigenvalues 



N 



d 



z.l„.\ 2^n z (L y ,d) = [ y 
d=o V h v 



(2.6) 



Since nz(L y , L y ) = 1, i.e., there is only a single \a,l v ,d,j 
for d = L y , we denote it simply as \a,l v ,l v - The single 
reduced eigenvalue with d = L y is 



Aa. 



L v ,Z/ 1 



= 1 



(2.7) 



We now proceed with our results. We shall point 
out relevant features for the widths that we consider; of 
course, it is possible to study larger widths, but, as our 
discussion will show, relevant features are already present 
for the widths that we consider. 



III. STRIPS OF THE SQUARE LATTICE 
A. Free 

We denote these strips as sq[L y x m],free. For L y — 
1, an elementary calculation yields Z = Q(Q + v) m ~ l . 
Setting Q = v 2 yields 



Z(sq[l x m],free,v z ,v) = v m+ \v + 1) 



(3.1) 



which has zeros only at the two discrete points v = 
and v = — 1. In this case the continuous B degenerates 
to these two points. 

For L y = 2 we use the partition function Z(sq[2 x 
m], free, Q, v), calculated in Ref. 01 1 which has the 
form (|2.1() with two A's. For Q — v 2 , we find 



A S( },2,0,j 



(v + 2) 2 ± (v 4 + Av 3 + 12v 2 + 20v + 12) 1/2 

(3.2) 



where j — 1,2 correspond to ±. In the infinite-length 
limit, the continuous accumulation set loci B v and, cor- 
respondingly, Bq, are shown in Figs. ^ an d [21 They 
consist of two complex-conjugate arcs that intersect each 
other and cross the real v axis at v = —2 and equiva- 
lently, the real Q axis at Q = 4. The endpoint of these 
arcs in the v plane occur at the roots of the polynomial 
in the square root of A S9j 2,0j : namely v — —1.33 ± 0.23« 
and -0.67 ± 2.48i, and hence Q = 1.71 ± 0.61i and 
Q = —5.71 ± 3.34i. For comparison, in this and other 
figures we show zeros of the partition function for long 
finite strips; in this case, m = 40. One sees that the zeros 
lie rather close to the asymptotic loci B and that the den- 
sity of zeros increases as one approaches the endpoints of 
the arcs. 

In these figures and others shown below, there are also 
zeros of the partition function that do not lie on the 
asymptotic accumulation loci. For example, in general, 
for any graph G, the cluster relation given above shows 
that Z(G,Q,v) = at the point (Q,v) = (0,0), which 
lies on manifolds defined by all of eqs. (|1.4|) - (|1.6|l . De- 
pending on the type of lattice strip graph, this may or 
may not be an isolated zero or lie on the continuous ac- 
cumulation set of zeros, B. For example, for the L y = 2 
square-lattice strips with free or cylindrical boundary 
conditions it is isolated (cf. Figs. Ill Hi . while for the 



L ,, 



2 strips with cyclic or toroidal boundary condi- 



tions, it lies on the loci B in the v and Q planes (cf. 
Figs. 151811 and similarly for the triangular strips to be dis- 
cussed below. Another general result is that for a graph 
G with at least one edge, Z(G,Q,v) = at the point 
(Q,v) = (1,-1). This follows because Z(G,Q,-1), the 
partition function for the zero-temperature Potts antifcr- 
romagnet, is precisely the chromatic polynomial P(G, Q), 
which counts the number of ways one can assign colors 
from a set of Q colors to the vertices of G, subject to the 
condition that no two adjacent vertices have the same 



3 



color. These are called proper colorings of G. Clearly, 
the number of these proper colorings of a graph G van- 
ishes if it has at least one edge and there is only one 
color, i.e., if Q = 1. This point (Q,v) = (1,-1) is on 
the manifold defined by eq. (|1.4I) for the square lattice 
(although not on the corresponding manifolds defined by 
eqs. Ijl.51) an( i <|1-6|> for the triangular and honeycomb 
lattices). Thus, one sees a zero at this point in the plots 
for the square-lattice strips. In the cases we have studied, 
this zero is isolated. 
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FIG. 1: Locus B v for the Potts model on a 2 X oo strip of the square 
lattice with free boundary conditions and with Q and v satisfying 
eq. 11.41 . Partition function zeros are shown for a 2 X 40 strip. 



For L y — 3 we use the exact calculation of Z(sq[3 x 
to], free, Q,v) for general Q and v in Ref. [l7T ] and spe- 
cialize to Q = v 2 . The partition function depends on 
(m'th powers of) four eigenvalues which are roots of a 
quartic equation (eq. (A. 8) in Ref. 17]). In the limit 
L x — ► oo, B v consists of complex-conjugate pairs of arcs, 
all of which pass through the point v = — 2, where all 
four roots of the above quartic equation are degenerate 
in magnitude. The arcs lying farthest from the real axis 
cross the imaginary v axis (at v ~ ±2.96i). Hence, the 
imagine of this locus under the map Q — v 2 in the Q 
plane, Bq also consist of complex-conjugate arcs which 
all pass through the point Q = 4. Furthermore, two 
of these cross the negative real Q axis, so that Bq sep- 
arates the Q plane into regions. There are also other 
zeros on the negative real v axis, and hence resultant ze- 
ros on the positive real Q axis. Us ing the calculations of 
Z(sq[L y x m], free, Q, v) in Ref. we have per- 

formed the corresponding analyses for L y — 4, 5 and have 
found similar features. 
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FIG. 2: Locus Bq for the Potts model on a 2 X oo strip of the 
square lattice with free boundary conditions and with Q and v 
satisfying eq. 11.41 . Partition function zeros are shown for a 2 X 40 
strip. 



B. Cylindrical 

We find that Z(sq[2 x m],cyl.,v 2 , v 
depending on two A's, which are 



has the form l|2.1() 



^sqcylj 



3v 2 + 8v + 6 ± (v + 2)\/bv 2 + 12^ + 8 



(3.3) 



where j — 1,2 correspond to ± and sqcyl refers to this 
type of strip. In the limit m — > oo, the locus B v consists 
of a self-conjugate arc that crosses the real axis at v = —2 
and has endpoints at the roots of the polynomial in the 
square root in eq. I|3.3|l . namely, v — (—6 ± 2i)/5. Thus, 
Bq consists of an arc that crosses the real axis at Q = 4 
and has endpoints at Q = (32 ± 24i)/25. These loci, 
together with partition function zeros, are shown in Figs. 
13 and 2| As was the case with the free strips, the density 
of zeros increases as one approaches the endpoints of the 
arcs. 



C. Cyclic and Mobius 

For Ly = 1, sq[l x to] is just the circuit graph 
with m vertices, C m - An elementary calculation yields 
Z(C m ,Q,v) = {Q + v) m +c^v m , so for Q = v 2 , one has 



A s <?,i.o =v + l and A sg ,i.i — 1 as in eq. I|2.7|l . and 
Z(C m ,v 2 ,v) =v m (v + l)[(v+l) m - 1 +v-l] . 
The resultant locus B v is the circle \v + 1| = 1, i.e. 



-1 



< d> < 2tt 



(3.4) 



(3.5) 



4 




for < <fi < 2tt. This locus crosses the real Q axis at Q — 
Q2 = 0, where it has a cusp, and at Q = Qi = Qoo = 4; 
it also crosses the imaginary Q axis at Q — ±2i. The loci 
B v and Bq divide the respective v and Q planes each into 
two regions. In the v plane these can be labelled as R\ 
and i?2, the exterior and interior of the circle \v + 1| = 1, 
and similarly in the Q plane the exterior and interior of 
the closed curve given by eq. Q3.6p . In regions R\ and 
i?2 the dominant A's are X S q,ifi and X S q,i,i, respectively. 

For L y — 2, we use the calculation of Z(sq[2 x 
m],cyc.,Q,v) in Ref. [0j- From cq. 1|2.4[1 we have 
n z (2,0) = 2 and n z (2,l) = 3, together with n z (2,2) 
1 . for a total of N Z ,2,X = 6. Specializing to the manifold 
of eq. I|1.4|l . we find that 



A 



sq,2,l,l 



= 1 



(3.7) 



A 



^9,2,1 ,j 



v + 2 ± V2v + 3 



(3.8) 



FIG. 3: Locus B v for the Potts model on a 2 X 00 strip of the 
square lattice with cylindrical boundary conditions and with Q 
and v satisfying the 2D phase transition condition 11.41 . Partition 
function zeros are shown for a cylindrical 2 X 40 strip. 



where j = 2, 3 correspond to ±, and X S q.2,o.j, given by 
eq. 
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FIG. 5: Locus B„ for the Potts model on a 2 X 00 cyclic or Mobius 
strip of the square lattice with Q and v satisfying eq. 11.41 . Parti- 
tion function zeros are shown for a cyclic 2 X 40 strip. 



FIG. 4: Locus Bq for the Potts model on a 2 X 00 strip of the 
square lattice with cylindrical boundary conditions and with Q 
and v satisfying eq. 11.41 . Partition function zeros are shown for a 
cylindrical 2 X 40 strip. 



which crosses the real v axis at v = and v = — 2. The 
resultant locus Bq with Q — v 2 is given by 

Re(Q) = 2cos0 (cos</> — 1) , 

Im{Q) = 2sin</) (cos0- 1) (3.6) 



For m — * 00, the locus B v for this cyclic (or corre- 
sponding Mobius) strip, shown in Fig. [S] is comprised 
of a single closed curve that intersects the real v axis at 
v = and v = —V2 and in a two- fold multiple point at 
v = —2. The image of B v in the q plane, Bq, shown in 
Fig. |H| is again a closed curve that crosses the real Q 
axis once at Q — and Q = 2 and in a two-fold multi- 
ple point at Q = 4, separating the complex Q plane into 
three regions in 1-1 correspondence with those in the v 
plane. These regions are 

• containing the real intervals v > and v < —2 
and extending outward to the circle at infinity, and 



5 



separate the respective v and Q planes in several regions 
in which different Ag^^^'s are dominant. 



10- 




y e- 

V lm(Q) 

4- 


V 






-8 -4 -a^"*' 


|w^4N> 6 8 


oj^^ -2- 


Re(Q) Y 


-4- 




X -6- 
° V "^--^_-8- 




-10- 





We have performed corresponding calculations for the 
cyclic and Mobius strips of the square lattice with L y = 4 
and L y — 5. For brevity, we only comment on Bq 
here. We find that Bq crosses the real axis in the in- 
terval 0<Q<4atQ = 4 and at Q = Q21 for integer 
1 < I < L y and also on the negative real axis. The 
outermost complex-conjugate curves on Bq continue the 
trend observed for smaller widths, of moving farther away 
from the origin. For example, the outermost curves on 
Bq cross the imaginary Q axis at Q — ±2i for L y = 1, 
Q ~ ±8.5i for L y = 2, and at progressively larger values 
for larger L y . Similarly, this outermost curve crosses the 
negative real axis farther away from the origin; the ap- 
proximate crossing point for L y — 3 is at Q ~ —10, with 
larger negative values for L y = 4, 5. 



D. Toroidal and Klein-bottle 



FIG. 6: Locus Bq for the Potts model on a 2 X 00 cyclic or Mobius 
strip of the square lattice with Q and v satisfying eq. 11.41 . Parti- 
tion function zeros are shown for a cyclic 2 X 40 strip. 



its image in the Q plane, containing the real inter- 
vals Q > 4 and Q < 0, in which (with appropriate 
choice of the branch cut for the square root in eq. 
(|3.2() ) A sgi 2,o,i is dominant, 

• i?2, containing the real interval — \f2 < v < 0, and 
its image in the Q plane containing the real interval 
< Q < 2, in which A S(?j 2,i,2 is dominant, 

• i?3, containing the real interval — 2 < v < — \[2 and 
its image in the Q plane, containing the real interval 
2 < Q < 4, in which X sq ,2,2 = 1 is dominant. 

As was the case for the cyclic L y — 1 strip, the curve Bq 
has a cusp at Q = 0. For comparison with the asymptotic 
loci, in Figs. [SJandElwe also show partition function zeros 
calculated for a long finite strip, with m = 40. One sees 
that these lie close to the respective loci B. 

The exact calculation of Z{sq[i x to], eye, Q, v) in Ref. 
[13 has the form of eq. (|2~3|) with L y = 3. From eq. (|2~4l) 
we have n^(3, 0) = 5 and 71^(3,1) = 9, 71^(3,2) = 5, 
together with 71^(3,3) = 1, for a total of Nz,z,\ = 20. 
We specialize to the manifold in eq. (|1.4|l . Owing to 
the large number of A sgj 3,dj's, we do not list them here. 
In the m — > 00 limit, we find that B v crosses the real 
v axis at v — v-i — 0, v — V4 = —y2, v = vq = — V3, 
and v = — 2, enclosing several regions in the v plane. The 
image locus under the map 1)1. 4f> . Bq, crosses the real axis 
in the interval < Q < 4 at Q = Q 2 = 0, Q = Q 4 = 2, 
Q = Qe = 3, and Q — 4. It also crosses the negative 
real axis at two points corresponding to the two pairs of 
complex-conjugate points away from v — at which B v 
crosses the imaginary axis in the v plane. As with the 
L y = 1 and L y = 2 cyclic strips, the curves on Bq and B v 



The exact solution for the partition function on the 
L y = 2 strip with toroidal boundary conditions which we 
obtained in Ref. 01 nas the form of eq. (|2.1|) with six 
A's. Setting Q = v 2 (and using the abbreviation sqtor to 
indicate the boundary conditions), we find 



A 



sqtor,2,j 



= A 



sqcyL2,j 



(3.9) 



where the \ S qcyi,2,j with j = 1,2 were given above in eq. 

El. 



A 



sqtor,2,j 



{v + 2)(v + 3) 



± 



(v 2 + 3v + 8 + 4V2 ) (v 2 + 3v + 8 - 4y/2 ) 



where j = 3, 4 correspond to the ± signs, 

A S( jtor,2,5 = t> + 1 , 

and 



A. 



= 1 



1/2' 

(3.10) 
(3.11) 
(3.12) 



In the m — > 00 limit of this strip with toroidal or Klein- 
bottle boundary conditions, we find that B v , shown in 
Fig. intersects the real v axis at v = v — —y2, and 
v = —2. The image curve Bq, shown in Fig. thus in- 
tersects the real axis at Q = 0, 2, 4. These curves divide 
the respective v and Q planes into three regions. In the 
v plane, these are (i) the region Ri including the semi- 
infinite intervals v > and v < —2 and extending to com- 
plex infinity, in which (with appropriate choices of branch 
cuts for the square roots) X sq t or ,2,i is dominant; (ii) the 
region i? 2 including the real interval — v2 < v < 0, en- 
closed by the outer curve, in which X sq tor,2,3 is dominant, 
and (iii) the region R3 enclosed by the innermost curve 



G 




FIG. 7: Locus B v for the Potts model on a 2 X oo strip of the square 
lattice with toroidal or Klein-bottle boundary conditions and with 
Q and v satisfying eq. 11.41 . Partition function zeros are shown for 
a toroidal 2 X 40 strip. 
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FIG. 8: Locus Bq for the Potts model on a 2 X oo strip of the square 
lattice with toroidal or Klein-bottle boundary conditions and with 
Q and v satisfying 11.41 . Partition function zeros are shown for a 
toroidal 2 X 40 strip. 



and including the real interval — 2 < v < —y/2, in which 
A sg tor,2,6 = 1 is dominant. Corresponding results hold in 
the Q plane. 

Using our results in Ref . 0, |2(j , we have also per- 
formed similar calculations for the L y = 3 strip of the 
square lattice with toroidal boundary conditions. We 
find that B v crosses the real axis at v = 0, v — —\/2 1 



v = —2, and v ~ —5.2 and contains complex-conjugate 
curves extending to complex infinity in the Re{v) < 
half-plane. Again, corresponding results hold in the Q 
plane. 



IV. STRIPS OF THE TRIANGULAR LATTICE 
A. General 

In order to investigate the lattice dependence of the loci 
B, we have also calculated these for infinite-length strips 
of the triangular lattice with various boundary conditions 
and with Q and v restricted to satisfy the phase transition 
condition for the two-dimensional triangular lattice, eq. 
1)1.5(1 . We construct a strip of the triangular lattice by 
starting with a strip of the square lattice and adding 
edges connecting the vertices in, say, the upper left to the 
lower right corners of each square to each other. Since 
we will find that the values Q21 for £ = 1, • • • , L y play 
an important role for the loci Bq for these strips with 
periodic longitudinal boundary conditions, just as they 
did for the corresponding square-lattice strips, we give a 
general solution of eq. 1)1.5)1 for the case where Q = Q r : 
the roots of this equation are —1 + 2 cos(2(r + 77)71") with 
r] = 1, —1, 0, i.e., in order of increasing v, 



vti{r) = — 1 + 2 cos 



Vt2(r) = -1 + 2 cos 



2(r + 1)tt 
37 



2(r- 1)tt 
37 



and 



1 + 2 cos 



(4.1) 



(4.2) 



(4.3) 



The relevant case here is r — 11 with 1 < I < L y . More 
generally, for any real r > 2, these roots have the prop- 
erties 



vn(r) < v t 2{r) < for r > 2 

(where the first equality holds only at r = 
second equality holds only at r = 2), 

v t s(r) > for r > 2 



(4.4) 



00 and the 



(4.5) 



(where the equality holds only at r = 2). As r increases 
from 2 to 00, (i) Vt\(r) increases from —3 to —2, (ii) Vt2{r) 
decreases from to —2, and (iii) Vta(r), which is the 
physical root for the phase transition on the triangular 
lattice, increases from to 1. 

For the cases of interest here, with r = 2£ and 1 < I < 
L y , for widths up to L y = 5, many of the trigonometric 
expressions in eqs. 1)4. 1)1 - 1)4.3)1 simplify considerably, to 
algebraic expressions, and in some cases to integers, so it 
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is worthwhile displaying these roots explicitly. For r 
and r — 4 we have 



i(2) = -3, va(2)=0, wt 3 (2)=0 =► 
-1 - v/3, v«(4) = -1, «t 3 (4) 



«ti(4) 



o = q 2 = o 

(4.6) 



-V3 
(4.7) 



-1 



We list the expressions for the Uy(r), j = 1,2,3 for r 
6, 8, 10 in the Appendix. For r = oo, 



«ti(oo) = w t2 (oo) 

Q — Qoo 4 



-2, u t3 (oo) = 1 



(4.8) 



(Outside of our range, at r = 1, since Qi = Qoo, eq. 
(|1.5|) with Q = Qi has the same set of roots as in 1)4. 8p . 
with uti(l) = w*3(l) = —2 and Wt2(l) = 1-) Concerning 
the behavior of the solutions of cq. (|1 . f>f> with v as the 
independent variable, as v decreases from 0, Q increases 
from 0, reaching a maximum of 4 at v = —2 and then 
decreasing through to negative values as v decreases 
through —3. Thus, for all v in the interval — oo < v < 0, 
Q is bounded above by the value 4. 
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FIG. 9: Locus B v for the Potts model on a 2 X oo strip of the 
triangular lattice with free boundary conditions and with Q and v 
satisfying cq. 11.51 . Partition function zeros are shown for a free 
2 X 20 strip. 



Free 



The exact solution for the partition function on the 
L y = 2 strip of the triangular lattice with free boundary 
conditions in Ref. [21| has the form of eq. <|2.1|1 with two 
A's. Setting Q — v 2 (v + 3) as in eq. I|1.5p . we find the 
corresponding reduced A's 



A 



t,2,0,j 



(V + 1) 



v 3 + 5v 2 + 9v + 7 



± (v + 3) ^{v + l)(i> 3 + 3v 2 + 3v + 5) 



(4.9) 



where j = 1, 2 correspond to the ±. 

In the limit m — > oo, we find the locus B v shown in Fig. 
consisting of the union of (i) a curve that crosses the 
real v axis at v = — 2 (a multiple point on the curve) and 
v = —3 and has cndpoints at the two complex-conjugate 
roots of the cubic factor in the square root in eq. 14.911 . 
v ~ —0.206 ± 1.37i and (ii) a line segment on the real 
v axis extending from v = — 1 to v ~ —2.587. These 
endpoints of the line segment are the other two zeros of 
the polynomial in the square root in eq. I|4.9|l . This locus 
divides the complex v plane into two regions. The image 
of this locus under the mapping (|1.5I) . Bq, is shown in 
Fig. 1 101 and consists of the union of a closed curve passing 
through Q = and Q = 4 with complex-conjugate arcs 
passing through Q — 4 and terminating at endpoints 
Q ~ —4.38 ± 4.12i, and a line segment on the real axis 
extending from Q = 2 to Q ~ 2.762. In Figs. Hand ESI 
we also show zeros calculated for a long finite free L y = 2 




FIG. 10: Locus Bq for the Potts model on a 2 X oo strip of the 
triangular lattice with free boundary conditions and with Q and v 
satisfying eq. 11.51 . Partition function zeros are shown for a free 
2 X 20 strip. 



strip of the triangular lattice. One sees again that these 
lie close to the asymptotic loci B v and Bq, as one would 
expect for a long strip. There are also discrete zeros 
that do not lie on (or close to) the equimodular curves 
B. We have already explained the origin of the zero at 
(Q,v) — (0,0). In contrast to the situation with free 
strips of the square-lattice, where the zero at v = and 
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its image at Q = were both isolated, here we find that 
for the free L y = 2 strip of the triangular-lattice (see Figs. 
l9l HOjl. there is an isolated zero at v = 0, but because of 
the non one-to-one nature of the mapping l|1.5l) , its image 
in the Q plane is not isolated but rather is on Bq. We 
have also performed analogous studies of wider strips of 
the triangular lattice using the calculations of Ref. (2^ j 
and we find qualitatively similar results. 

C. Cyclic and Mobius 

For cyclic boundary conditions we have Nz,2,\ = 6, 
n z (2,0) = 2, n z (2,l) = 3, n z (2,2) = 1, just as in the 
square-lattice case. We calculated the general partition 
function Z(t[2 x to], eye, Q, v) in Ref. [2J|, and this has 
the form of eq. (|2.3|) with L y _ — 2. Restricting Q and v to 
satisfy eq. i|1.5|) . we obtain \t,2,2 = 1, M,2,o,j f° r J = 1,2 
given in eq. (|4.9(l . and Xt,2,i,j for j = 1,2,3, which are 
solutions to the equation 

jf - (v + 2)(3v + 4)?7 2 + (v + 2)(v + l)(3v 2 + 9v + 4)rj 



-{v+l) 2 (v 2 + 3v + lf = . 


(4.10) 




r2 




-1 im(v) 

; 




, , , , 




—1 




—2 



FIG. 11: Locus B v for the Potts model on a 2 X oo strip of the 
triangular lattice with cyclic or Mobius boundary conditions and 
with Q and v satisfying eq. 11.51 . Partition function zeros are 
shown for a cyclic 2 X 20 strip. 

In the limit to — > oo for this cyclic strip (and for the 
same strip with Mobius boundary conditions), we find 
the locus B v shown in Fig. ^] and the corresponding 
locus Bq shown in Fig. ^] These loci consists of closed 
curves that intersect the real v and Q axis at the following 
points, where the value of Q is the image of the value of 
v under the mapping (|1.5|l : (i) v = and v — —3, =>■ 
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-10 



FIG. 12: Locus Bq for the Potts model on a 2 X oo strip of the 
triangular lattice with cyclic or Mobius boundary conditions and 
with Q and v satisfying eq. 11.51 . Partition function zeros are 
shown for a cyclic 2 X 20 strip. 

Q = Q2 = 0, (ii) v = -1, => Q = Qi = 2, and (iii) 
v = —2, =>• Q = 4. The locus B v also contains a line 
segment extending from v ~ —1.531 to v = —2, and Bq 
contains its image under the mapping Q1.5[l . extending 
from Q ~ 3.44 to Q = 4. 

The locus separates the i> plane into four regions, 
as is evident in Fig. 1111 

• the region including the intervals v > and v < —3 
on the real v axis and extending to complex infinity, 
in which (with appropriate definition of the branch 
cuts associated with the square root in eq. I|4.9|)) 
At, 2, 0,1 is dominant, 

• the region including the neighborhood to the left 
of v = and excluding the interior of the loop 
centered approximately around v = —1.2, in which 
one of the At,2,i.j's is dominant, 

• the region in the interior of the loop centered 
around v — —1.2, in which Xt,2,2 = 1 is dominant, 

• the region including the interval that extends from 
v = —2 to v = —3, in which another \t,2,i,j is 
dominant. 

We give some specific dominant eigenvalues at special 
points. At v = 0, ^{,2,0,1 = (7 + 3v / 5)/2, equal in mag- 
nitude to the dominant At,2,i,j- At v = —2, all of the 
six \t,2,d,j's are equal in magnitude, and equal to unity. 
At v = —3, I At 2,o,j 1 — 2 for j — 1,2, and two of the 
I A* 2 i,j I = 2, while the third is equal in magnitude to 1. 

Correspondingly, the locus Bq divides the Q plane into 
five regions: 
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• the region including the intervals Q > 4 and Q < 
on the real Q axis and extending to infinity, in 
which At,2,o,i is dominant 

• two complex-conjugate regions bounded at large Q 
by curves that cross the imaginary axis at Q ~ 
±8.23*, in which one of the At,2,i,j's is dominant 

• the region including the interval < Q < 2, in 
which another Xt.2,i.j is dominant 

• the region in the interior of the loop centered ap- 
proximately around Q = 2.6, in which At : 2,2 = 1 is 
dominant 

We have performed similar calculations for cyclic strips 
of the triangular lattice with greater widths, L y = 3, 4, 5. 
We find that B v crosses the negative real axis at «ta(2^) 
for 1 < £ < L y , so that Bq crosses the real axis at the 
image points under eq. H1.5J) . Q2i- As with the cyclic 
square-lattice strips with widths L y > 3, we find that the 
outermost curves on Bq cross the negative real axis; for 
example, for L y = 3, such a crossing occurs at Q ~ —9.4. 

D. Other Strips of the Triangular Lattice 



and (iii) Vh c s(r), the physical root, increases monotoni- 
cally from to 4. As with the analogous expressions for 
the triangular lattice, it is straightforward to work out 
simpler expressions to which eqs. (|5.1[1 - H5.3[) reduce for 
special values of r; we omit the details here. 

For this honeycomb lattice, substituting the value v = 
—2 into eq. I|1.6fl yields two solutions, Q = 2 and Q = 4. 
Using our general solutions for the partition functions 
Z(hc\2 x m], free,Q,v) and Z(hc[2 x m],cyc.,Q,v) in 
Ref. |15l| . we have checked that at v = — 2, Q = 4 there 
is degeneracy of dominant A's, so that these points are 
on the respective loci B v and Bq. This property is thus 
similar to the feature that we have found for strips of 
the square and triangular lattices. For L y — 2 cyclic 
honeycomb-lattice strips we find that, in addition, the 
points v — 0, Q = Qi — and v — —2, Q — Q4 — 2 
are on the respective loci B v and Bq. For the L y — 3 
cyclic strips of this lattice, B v contains these points and 
also Vhd(6) = — 2\/3 cos(57r/18), so that the image Bq 
contains the point Q = Qq — 3. For all the strips of 
this lattice that we have studied, with L y up to L y = 5, 
we find that B v crosses the real axis at v = vi lc i(2£) for 
£ = 1, • • • , L y , so that Bq crosses the real Q axis at Q21 
for £ = !,••• , L y . 



We have also calculated the partition function and the 
resultant loci B v and Bq for strips of the triangular lat- 
tice with cylindrical and toroidal boundary conditions. 
A general feature that we find is that B v passes through 
v = —2, and hence Bq passes through Q = 4. Other 
features depend on the specific boundary conditions and 
width. One property that we encounter is noncompact- 
ness of B v and Bq (as was the case with the L y = 3 
toroidal strip of the square lattice and for cyclic self-dual 
strips of the square lattice Q). 



V. STRIPS OF THE HONEYCOMB LATTICE 

We first give a general solution for the three roots in 
v of eq. with Q = Q r ; in order of increasing value, 

these are 



Vhd(r) = -4cos(^ cos 
V r 



7T /l 

3 Vr 



1 



Vhc2(r) 



-4 cos ( 





1 






j COS 








3 


vr 





and 



Vhdi(r) = 4cos( '— J cos(^ 



V r 



(5.1) 



(5.2) 



(5.3) 



As r increases from 2 to 00, and Q r thus increases from 
to 4, (i) Vhci{r) decreases from to a minimum of —9/4 
at r = 7r/ arcsin(-\/T0/8) ~ 7.73 and then increases to 
—2; (ii) vi l c2{r) decreases monotonically from to —2; 



VI. DISCUSSION 

We remark on several features that are common to all 
of the strips of the three types of lattices that we have 
analyzed. These include the following: 

• For all of the strips, including those of the square, 
triangular, and honeycomb lattices, that we have 
studied where nontrivial continuous accumulation 
loci are defined (thus excluding the lxm line graph 
with free boundary conditions), we find that B v 
passes through v = — 2 and Bq passes through Q = 
4, which is the image of v = —2 under both of the 
mappings (|1.4I) and l|1.5fl and is a solution of eq. 
!p~B|) with v = -2. 

• For all of the cyclic (and equivalently, Mobius) 
strips that we have studied, besides the crossing 
at Q = 4, Bq crosses the real axis at 



!2C 



for £ = 1, 



(6.1) 



We conjecture that this holds for arbitrarily large 
L y . This locus can also cross the real axis at other 
points, such as the crossings on the negative real 
axis that we found for widths L y > 3. 

For the particular case of the cyclic square-lattice strips, 
these properties agree with a result in Ref. |l3| , namely, 
that at Q = Q2e, \X S q,L y j,max\ = |A sg: L y /_i : „ja2;|, where 
K q ,L v ,i,raax denotes the eigenvalue X sqj L y ,£,j of largest 
magnitude. 

We also observed that the points v = — 2 and Q = 
4 play a special role for self-dual strips of the square 
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lattice in Ref. [4j . It is interesting that for self-dual cyclic 
square-lattice strips, in addition to the point Q = 4, Bq 
crosses the real axis at Q^+i f° r I < ( < L y . This 
set of points is interleaved with those in eq. 1)6.10 . As 
we have noted in Ref. [4|, these findings are in accord 
with the fact that the Potts model at the values Q — 
Q r has special properties, such as the feature that the 
Tempe rley -Lieb algebra is reducible at these values 0, 
EJ El E| ■ In Ref- we compared our exact results for 
B on infinite-length, finite-width cyclic self-dual strips of 
the square lattice with calculations of partition function 
zeros for finite L x L sections of the square lattice with 
Q = v 2 in Ref. with the same boundary conditions. 
For finite lattice sections, there is, of course, no locus 
Bq defined, and hence one is only able to make a rough 
comparison of patterns of zeros. In the calculation of 
zeros for the above-mentioned Lx L section of the square 
lattice with cyclic self-dual boundary conditions, e.g., for 
L = 8, a number of zeros in the Q plane occur at or 
near to certain Q^s, and the zeros in the v and Q planes 
exhibit patterns suggesting the importance of the points 
V = —2 and Q = 4. With our results in the present 
paper, we can extend this comparison. We see that the 
importance of v = —2 and Q — 4 for the pattern of 
partition function zeros for Q and v satisfying the relation 
(|1.4|) generalizes to square-lattice strips with a variety 
of boundary conditions, not necessary self-dual. Indeed, 
going further, our results show that the features we have 
observed are true not just of the partition function on 
the square-lattice strips with Q = v 2 , but also on strips 
of the triangular and honeycomb lattices with Q and v 
satisfying the analogous phase transition conditions (|1.5|) 
and (|l.(j|) . One interesting aspect of the findings in the 
present work and our Ref. [4| on infinite- length, finite- 
width strips is the special role of the value Q — 4 for 
the loci Bq, which can make a connection with the locus 
B for the physical phase transition of the Potts model 
on two-dimensional lattices. In this context, we recall 
that the value Q — 4 is the upper end of the interval 
< Q < 4 for which the ferromagnetic Potts model has 
a second-order transition on two-dimensional lattices. 



VII. CONCLUSIONS 

In conclusion, we have presented exact results for the 
continuous accumulation set B of the locus of zeros of 
the Potts model partition function for the infinite-length 
limits of strips of the square, triangular, and honeycomb 
lattices with various widths, a variety of boundary con- 
ditions, and Q and v restricted to satisfy the conditions 
(|1.4|) , 1)1. 5|) , and l|1.6f) for the ferromagnetic phase transi- 
tion on the corresponding two-dimensional lattices. We 
have discussed some interesting general features of these 
loci. 
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IX. APPENDIX 

In this Appendix we list the Vtj{r) for r = 6, 8, and 

10, which are relevant to the analysis of the cyclic strips 

of the triangular lattice for widths up to L y — 5. From 
eqs. I|4.1|) - (I4.3|) . we have 

«y(6), J = 1,2,3 Q = Q 6 -3 (9.1) 

where [33| 

uti(6) = -1 + 2cos(7tt/9) ~ -2.532089 (9.2) 

ut2(6) = -1 + 2cos(5tt/9) ~ -1.347296 (9.3) 

w*3(6) = -1 + 2cos(tt/9) ~ 0.879385 (9.4) 

v tj (8),j = 1,2,3 =^ Q = Q 8 = 2 + V2 (9.5) 
where 

v t i(S) = -l-V2~ -2.414214 (9.6) 



1 — \/3 

«ta(8) = -l+ ^---1.517638 (9.7) 

v2 



M8) = -1 + X - 0.931852 
v2 



u ti (10), j = 1,2,3 ==► Q = Q 10 = ^±^L (9.9) 
where 

uti(10) = ^ (~ 5 + V5 - y30 + 6%/5 J ~ -2.338261 



(9.10) 



(1 + V5) ^ L 



v t2 (W) = ' v ~ ' ~ -1.618034 (9.11) 



uj 3 (10) = i (-5 + a/5 + y30 + 6\/5 j ~ 0.956295 



(9.12) 
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